Metamagnetic phase transition of the antiferromagnetic Heisenberg icosahedron 
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The observation of hysteresis effects in single molecule magnets like Mni2-acetate has initiated 
ideas of future applications in storage technology. The appearance of a hysteresis loop in such 
compounds is an outcome of their magnetic anisotropy. In this Letter we report that magnetic 
hysteresis occurs in a spin system without any anisotropy, specifically, where spins mounted on the 
vertices of an icosahedron are coupled by antiferromagnetic isotropic nearest-neighbor Heisenberg 
interaction giving rise to geometric frustration. At T = this system undergoes a first order 
metamagnetic phase transition at a critical field between two distinct families of ground state 
configurations. The metastable phase of the system is characterized by a temperature and field 
dependent survival probability distribution. 

PACS numbers: 75.50.Xx;75.10.Hk;75.10. Jm;75.30.Cr 
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Introduction — Low-dimensional magnetic systems 
show a variety of fascinating phenomena that are 
associated with geometrical frustration P, Among 
them are magnetization plateaus and jumps as well as 
unusual susceptibility minima, as observed for example 
for the kagome lattice antiferromagnet 0, 0| . Some of 
these effects can also occur in certain strongly frustrated 
magnetic molecules such as the Keplerate {Mo72Fe3o}|3- 
In this Letter we report that a first order metamag- 
netic phase transition (with associated hysteresis and 
metastability effects) occurs for a system of spins that 
are mounted on the vertices of an icosahedron when 
an external magnetic field, B, equals a critical value 
Be- These spins interact with one another only via 
nearest-neighbor, antiferromagnetic isotropic Heisen- 
berg exchange, but due to the geometrical frustration 
originating in the unique geometry of the icosahedron 
they undergo the metamagnetic transition despite the 
absence of any anisotropic energy terms. As the field 
proceeds through a closed cycle the magnetization traces 
the hysteresis loop shown in Fig. ^ The metastable 
phases, including the temperature dependence of their 
lifetime, show rich characteristics. 

Our exact classical treatment shows that the abrupt 
transition at T = originates in the intersection of two 
energy curves belonging to different families of spin con- 
figurations that are ground states below and above the 
critical field. The minimum of the two energy functions 
constitutes a non-convex minimal energy function of the 
spin system and this gives rise to a metamagnetic phase 
transition Q. At T = the partition function is a non- 
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FIG. 1: Hysteresis behavior of the classical icosahedron in 
an applied magnetic field obtained by classical spin dynamics 
simulations. 

analytic function of B, and since the magnetization fea- 
tures a finite jump at B^ the transition is of first order. 
We also show that the corresponding quantum spin sys- 
tem for sufficiently large spin quantum number s pos- 
sesses a non-convex set of lowest energy levels when plot- 
ted versus total spin. This is the discrete analog of the 
non-convex classical minimal energy function. Therefore, 
also the quantum spin system features an unusual mag- 
netization jump. 

Classical simulations at T ~ — The behavior of clas- 
sical spin systems subject to an applied magnetic field 
both at T « and finite temperatures can be very effec- 
tively studied with the help of a stochastic spin dynam- 
ics approach such as that proposed in 0. Here, the spin 
system is coupled to a heat bath in a Langevin-type ap- 
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proach by using a Landau-Lifshitz-damping term as well 
as a fluctuating force with white noise characteristics. 
Starting from an arbitrary initial configuration the spin 
system can be investigated either at zero temperature by 
observing the relaxation to its ground state or at finite 
temperature by following its time evolution. 

We consider first T = 0. For i? = we find, that the 
spins adopt a non-coplanar configuration, where each of 
the 12 spin vectors makes an angle of arccos(— 1/a/5) ~ 
116.6° with respect to its five nearest neighbors. This is 
in agreement with the analytical results derived in 8] . 

The spins are also subjected to an external magnetic 
field that increases from zero linearly with time but at a 
very slow rate. The dynamical evolution of the spins is 
monitored as a function of time and the results are stored 
within an animation file [isj enabling direct visualization 
of the spin vectors. In the first stage the configuration of 
the spin vectors evolves continuously with B until sud- 
denly there is an abrupt change in their orientations. In 
the new configuration two spin vectors at opposite ends 
of a diameter of the icosahedron are aligned parallel to 
the field while all of the remaining 10 spin vectors sub- 
tend a common polar angle with respect to B and their 
azimuthal angles are uniformly spaced. With increasing 
B the common polar angle decreases monotonically while 
the azimuthal angles remain fixed until at saturation all 
spin vectors are parallel to B. 

Analytical results — For a classical spin system the 
ground states are defined as states s = (s(l), . . . , s{N)) 
which minimize the energy 

His,B) = Hois)-BY,Sz{^l) (1) 

The coupling J^j^ between spins is chosen as J = 1 
between nearest-neighbors and zero otherwise. Accord- 
ingly the magnetic field is given in appropriate units. 

The (degenerate) ground states must fulfill the follow- 
ing necessary condition, compare Eq. (21) in 

J2J^Mi^)^^M + Ib = l,...,iV) , (2) 

V 

where denote suitable Lagrange parameters. Al- 
though this system of equations can only be solved nu- 
merically in most cases, we are confident that the follow- 
ing statements about ground states of the icosahedron 
subject to magnetic fields are correct 0. 

As remarked above, the spin configuration for B — 
has been established [Sj . The orientation of the individual 
spins is such that they form four groups of three spins 
where each group i is characterized by a common polar 
angle 9i and uniformly spaced azimuthal angles. 

For B > we numerically solve ^ with the assump- 
tion that the azimuthal angles remain fixed and only the 



four polar angles vary [l^ . Thus we obtain a 1-parameter 
family (the "4-6'-family") of possible ground states. Inter- 
estingly, it provides a local minimum of the energy only 
for < M < 5.61441. One might assume that the 4-6*- 
family also provides a global minimum of energy, i.e. a 
ground state, for the same interval < M < 5.61441, 
but this is wrong: There exists a different 1-parameter 
family of solutions of which has a lower energy than 
the 4-6'-family for M > Ma ~ 4.92949. As remarked 
above, this family can be characterized by two spin vec- 
tors which are aligned parallel to B and 10 spin vectors 
with a common polar angle 9 and uniformly spaced rel- 
ative azimuthal angles. The end points of the latter 10 
spin vectors form a regular decagon, and we will call this 
set of states the "decagon family" . As discussed below, 
the decagon family provides a local minimum of the en- 
ergy for 0.54102 < M < 12. These families give rise to 
two convex curves in the E vs. M diagram, Ei{M) for 
the 4-^?- family and E2{M) for the decagon- family, which 
intersect at the point with coordinates {Mq,Eq), where 
Mo ~ 4.92949 and Eq ~ -11.150. 
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FIG. 2: Schematic representation of two minimal energy 
curves Ei[M) and E2[M), whose minimum -Emin(Af) is not 
convex. Their common tangent (solid line) has a slope of Be- 

We have strong numerical evidence that the mini- 
mum of the two curves provides the absolute minimum 
E^in{M) of Ho{s) for given M 0. The latter function 
Einin{M) is therefore not convex and this translates to a 
jump in the magnetization, AM = M2 — Mi, at a critical 
field Be, compare 0. One can identify B^ as the slope 
of the common tangent of the curves Ei{M) and E2{M) 
illustrated in Fig. El This construction is fully equiva- 
lent to the statement that the total energies of the two 
phases, Si (Mi) - BMi and £^2(^2) - BM2, are equal 
for B = Be- According to the Ehrenfest classification the 
phase transition is of first order. Equivalently, B^ can be 
obtained by a Maxwell construction in the M vs. B dia- 
gram. The pertinent quantities for the phase transition 
are given in Tabled 

Since M{B) has a jump at B = B^ the susceptibility 
X = ^ diverges at B = Be and T = 0. It is rather easy 
to derive the leading temperature dependence of this di- 
vergence for T — !■ 0. Recall that an alternative expression 
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TABLE I: Characteristic values of the first-order phase transi- 
tion of the spin icosahedron, compare Fig.|21 -Bsat ~ 2(5 + \/5) 
denotes the saturation field. Xi = ^j^|s=Sc, i — 1,2 denote 
the limit values of the susceptibility. 



Be ~ 5.87614 
Ml ~ 4.71461 
Ei{Mi) ~ -12.4324 
Eo ~ -11.150 
Xi ~ 1.13294 



for the susceptibility is 



Be 0.40603 B^at 
Ma ~ 5.10784 
EiiMi) ~ -10.1218 
Mo ~ 4.92949 



X2 = 



0.8018 



X{B, P)^(3 {{M')p - (M)2) EE f3a^M) , (3) 

where ( )p denotes the thermal expectation value 
w. r. t. the canonical ensemble p = p{j3,B). In the 
case of a metamagnetic phase transition of the kind de- 
scribed above the ground state is degenerate for T ^ 
and B ^ Be- Hence a^{M) remains finite in this limit 
and x{Bc,P) diverges linearly with f3, i.e. the critical 
exponent is 1. 

Stability and hysteresis — In order to investigate the 
stability of the two families of ground states we performed 
a standard stability analysis by constructing the Hesse 
matrix, which for this system is a symmetric 24 x 24 ma- 
trix. Since the energy is invariant under rotations about 
the z-axis, one eigenvalue must be zero. We call a state 
satisfying Q "stable" if the Hesse matrix has only pos- 
itive eigenvalues apart from one zero eigenvalue. This 
is equivalent to the fact that the given state provides a 
local minimum of the energy. 

By applying this procedure to the two families of possi- 
ble ground states of the icosahedron we determine numer- 
ically the above-mentioned stability ranges of the fami- 
lies: The 4-6'-family is stable for < M < 5.61441 and 
the decagon family is stable for 0.54102 < M < 12. 

This implies that the system at T ~ will not imme- 
diately jump from the 4-0-family into the decagon family 
if B increases beyond B^ but remain in its family until 
M > 5.61441. Conversely, the decagon family will re- 
main the de facto spin configuration of the icosahedron 
if B is lowered beyond B^ until M < 0.54102. In fact, 
these hysteresis effects are observed in our simulational 
studies, see Fig. ^ 

Quantum calculations — We now discuss how the meta- 
magnetic phase transition is manifested in the quantum 
Heisenberg icosahedron. Classically, the phase transi- 
tion consists of a discontinuity of the magnetization as 
a function of the magnetic field. Quantum mechanically, 
the magnetization curve for T = is already a staircase 
of successive steps of unit height AM = 1, which re- 
sult from crossings of levels with adjacent total magnetic 
quantum numbers M and M -t- 1. In the context of this 
phase transition we are looking for a magnetization jump 



of unusual height, i.e. AM > 1. It is clear that such a 
jump must occur because it occurs in the classical limit 
s oo. The remaining questions therefore are, for which 
intrinsic spin quantum number s is such a jump clearly 
visible, and are there other signs of the phase transition 
at smaller s? 
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FIG. 3: Magnetization curves for T — and various values 
of the intrinsic spin quantum number s: The magnetization 
plateaus of smallest width are highlighted on each curve. At 
s = 4 a magnetization jump of AAI = 2 occurs, marked by the 
arrow. The solid curve shows that the field values that bisect 
the smallest plateaus converge to the classical transition field 
(dashed line). For large s the magnetization curve could only 
be evaluated down to a certain field due to prohibitively large 
dimensions of the related Hilbert spaces. 

Using a Lanczos procedure which yields numerically 
exact lowest energy eigenvalues in subspaces with con- 
stant total magnetic quantum number AI we are able 
to evaluate magnetization curves at T = for various 
intrinsic spin quantum numbers. Figure |31 shows the rel- 
evant parts of the magnetization curves at T = for 
s — 1/2, . . . , 4. For integer values of s the magnetization 
plateau of smallest width is highlighted; for half-integer 
values of s we highlight two such plateaus. With increas- 
ing s these widths shrink and already at s = 4 a magne- 
tization jump of AM = 2 occurs. This corresponds to a 
non-convex part of the discrete energy levels versus M. 

In Fig. |51 we also provide a curve which bisects the 
magnetization plateaus of smallest width. Assuming that 
the bisector value of B /i?sat is described by a polynomial 
in 1/s, we obtain as an estimate for the classic transition 
field B/Bc « 0.40±G.01, which is in very good agreement 
with the classical result (see Table . The uncertainty 
originates from the limited number of data points (eight) 
as well as from fluctuations between integer and half- 
integer values of s. 

Classical finite-temperature simulations — We have re- 
stricted our investigation to the determination of the 
lifetime of the high-field phase (decagon-family) in its 
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metastable regime {B/Bc < 1). The lifetime has been 
determined by the foUowing procedure: First, the sys- 
tem is prepared in the decagon phase at T = 0. Then 
the field is lowered to a value below the critical field value 
Be- Starting from these initial conditions the tempera- 
ture is set to a value T > and the trajectory of the 
system is calculated numerically by solving the stochas- 
tic Landau-Lifshitz equation. 

Shown in the inset of Fig. ^ is a single trajectory of 
a sample system. We exploit the unique property of the 
decagon family that two spins persist in pointing in the 
direction of B {sz{l) ~ Szi^) w 1) until abruptly break- 
ing away. This allows one to obtain an accurate determi- 
nation of the decay time for this system. By performing 
10^ such runs for each choice of T and analyzing the his- 
tograms of the resulting decay times one can determine 
the lifetime distribution. A common measure is the so- 
called survival probability Ps{t) which is the probability 
of the metastable state not having decayed by the time 
t. In Fig. ^ we have plotted Psit) for the metastable 
state for various temperatures and an external field in 
the metastable regime. An appropriate choice for the 
lifetime, tg, is the root of Ps{ts) = 0.5. We find that 
ts increases with decreasing temperature and appears to 
diverge for T — > as 1/T. Although one obtains sim- 
ilar probability distributions for systems showing ther- 
mally activated magnetization switching i we empha- 
size that our model Hamiltonian does not contain any 
additional energy term providing an energy barrier. In 
fact, it is the special geometry of the icosahedron that 
causes the system to show metastability. 
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FIG. 4: Survival probability Ps{t) for the metastable decagon 
phase subject to an external field B/Bsat = 0.27 for tem- 
peratures ksT/J = 0.025,0.015,0.005,0.0005 (left to right). 
Inset: Example trajectory of the system's time evolution at 
finite temperature. 

Summary and outlook — In this Letter we have shown 
that the antiferromagnetic Heisenberg spin icosahedron 
undergoes a metamagnetic phase transition and displays 
rich hysteresis and metastability phenomena when sub- 



ject to a varying external field. Given this, the metam- 
agnetic transition of the Heisenberg icosahedron may be 
of interest for potential applications in the area of nano- 
magnetic switches. It is also relevant in connection with 
magnetocalorics, since the magnetization jump is accom- 
panied by an enhanced magnetocaloric effect 0, ■ It 
is therefore very encouraging that recent progress in the 
synthesis of magnetic molecules offers t he p rospect of re- 
alizing the Heisenberg icosahedron 0, Q| . 
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